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Ge 2 Sb 2 Tes (GST) has been widely used as a popular phase change material. In this study, we 
show that it exhibits high Seebeck coefficients 200 - 300 pV/K in its cubic crystalline phase (c-GST) 
at remarkably high p-type doping levels of 1x10^® - 6x10^® cm~^ at room temperature. More 
importantly, at low temperature (T = 200 K), the Seebeck coefficient was found to exceed 200 pV/K 
for a doping range 1x10^® - 3.5x10^® cm“^. Given that the lattice thermal conductivity in this 
phase has already been measured to be extremely low (~ 0.7 W/m-K at 300 K),— our results suggest 
the possibility of using c-GST as a low-temperature thermoelectric material. 

PACS numbers: Computational Physics, Condensed Matter Physics 


I. INTRODUCTION 

Low-temperature thermoelectric (TE) materials are 
important for refrigeration applications in electronics, in¬ 
frared detectors, computers and other areasHowever, 
in contrast to high temperature TE materials, there is 
relatively little progress in the discovery of high perfor¬ 
mance low temperature TE materials. The most widely 
used materials at ambient temperature and below are 
(Bi,Sb) 2 Te 3 derived compounds^i^ These have figures of 
merit, ZT, of roughly unity at 300 K with lower ZT val¬ 
ues at lower temperatures. Here we show that cubic GST 
phase change material has potential to become an excel¬ 
lent thermoelectric material in the important tempera¬ 
ture range from 200 K - 300 K. Specifically, we find using 
Boltzmann transport calculations that disordered cubic 
p-type GST will have high thermopowers in the range 
200 /rV/K - 300 /K at remarkably high doping levels 
consistent with good conductivity. GST already is known 
to have a very low thermal conductivityii 

Pseudobinary alloys of (GeTe)m-(Sb 2 Te 3 )„ have been 
used in data storage for more than a decade because 
of their fast phase switching between metastable crys¬ 
talline (cubic) and amorphous phases A— The above 
room temperature thermoelectric properties of (GeTe)m- 
(Sb 2 Te 3 )„ were previously studied experimentallyAJ^ In 
particular, for m = 2 and n = 1; ordered Ge 2 Sb 2 Te 5 was 
reported to show interesting thermoelectric properties, 
even at room temperature ^ d ^ 

Disorder in the crystal structure can be beneficial for 
thermoelectric propertie s^^d'^ and hexagonal GST has 
been already reported to have reasonably good ther¬ 
moelectric properties4d^ The disorder in phase change 
materials was discussing in detail by Zhang and co- 
workersj^ We therefore investigate the thermoelectric 
properties of cubic disordered GST. Previous measure¬ 
ments of thermoelectric properties in c-GST were re¬ 


stricted to thin-film samples where the thermoelectric 
properties were found to vary significantly depending on 
the thickness of the sample liidi The thermoelectric prop¬ 
erties of c-GST in the bulk form have remained unknown. 
We note that Lee and co-workers have discussed resonant 
bonding in the context of thermoelectric performance^^ 
and that resonant bonding has been identified in phase 
change materials 


II. THEORETICAL APPROACH AND 
COMPUTATIONAL DETAILS 

The ab-initio calculations presented here were carried 
out using the linearized augmented plane-wave (LAPW) 
method^i as implemented in the WIEN2K code<^ The 
structures were relaxed using the PBE-GGA functional 
then a semilocal functional of Tran and Blaha (TB- 
m BJ}^ was employed to treat the exchange-correlation 
potential in the calculations of electronic, optical and 
transport properties. This functional form was shown 
to better describe the band gaps of a variety of semicon¬ 
ductors and insulatorsi ^^'^^ 

The cubic phase of GST occurs in a rock-salt struc¬ 
ture with Te occupying one site and Ge, Sb and vacancies 
randomly occupying the other site. We generated special 
quasirandom structures (SQS) using the prescription of 
Zunger et al^ for our first principles study. The SQSs 
were generated such that all the pair correlation functions 
(PGF) up to the next nearest neighbors are identical to 
the average PCF of an infinite random structure. We 
used the mcsqs code as implemented in ATAT software 
package^l to generate the SQSs using a Monte Garlo algo¬ 
rithm. We considered two SQS structures with 27 atoms 
(6 Ge, 6 Sb, and 15 Te atoms) and 45 atoms (10 Ge, 
10 Sb, and 25 Te atoms) in the unit cell and compared 
their electronic and transport properties. Despite the rel- 
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atively poor quality of PCFs for the case with 27-atomic 
unit cell, we find that the overall calculated properties for 
these structures are very similar (a detailed comparison 
of electronic, optical and transport properties in these 
structures can be found in Refi^). This may indicate 
that the inherent properties in this phase is independent 
of its PCFs. Thus this 45-atomic unit cell is converged 
with respect to size to calculate properties of c-GST. 

The self-consistent calculations were performed with a 
4x4x4 k-point grid of 36 k points in the irreducible Bril- 
louin zone (IBZ) whereas electronic and optical proper¬ 
ties were calculated using a 8x8x8 k-mesh. For the cal¬ 
culation of the thermoelectric properties, we employed 
a denser k-point mesh of 14x14x14. All the calcula¬ 
tions were duly tested for convergence with a variety of 
k-point grids. The thermopower was analysed using the 
Boltztrap code^S employing Boltzmann transport the¬ 
ory within the constant scattering time approximation 
(CSTA). This CSTA approach has been shown to be suc¬ 
cessful in calculating the Seebeck coefficient in a variety 
of thermoelectric materialsi^^— The substance of CSTA 
is to assume that the energy dependence of the scatter¬ 
ing rate is small compared with the energy dependence 
of the electronic structure. A detailed description of this 
approach can be found elsewhereWe used well con¬ 
verged basis sets with LAPW basis size corresponding 
to RMT^max=d-0, where Rmt is the smallest Muffin 
Tin radius and K^ax is the plane-wave cutoff parame¬ 
ter. The LAPW sphere radius of 2.47 Bohr was used for 
Ge whereas 2.5 Bohr was used for Sb and Te. Spin-orbit 
coupling was included for all the properties reported. 


III. RESULTS AND DISCUSSIONS 

A. Electronic structure and optical properties 

The total density of states (DOS) of c-GST near the 
Fermi energy (Ei?) is shown in Fig. [T] We find a clear 
semiconducting behavior with a band gap of 0.17 eV. 
This is in good agreement with a previous theoretical 
investigation^ where an indirect band gap of 0.1 eV was 
noted for c-GST although they used different exchange- 
correlation potential (local density approximation with¬ 
out spin-orbit coupling) and different method (shifting of 
the hexagonal structure) to generate the cubic structure. 
Another previous study considered the cubic GST as a 
narrow-gap degenerate semiconductor, with reside in¬ 
side the valence band or a defect bandit The experi¬ 
mentally reported estimated optical gap is 0.5 eV^i^s 
which is larger than what we obtain. 

GST has applications in rewritable DVD and Blu- 
ray technologies where a large contrast between the op¬ 
tical properties of the cubic and amorphous phases is 
exploitedj^ii In Fig. [2] and [3l we present the averaged 
reflectivity and optical absorption coefficients for c-GST, 
respectively. It is evident from Fig. [5] that c-GST dis¬ 
plays a high reflectivity ranging from 0.65 to 0.68 within 
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FIG. 1: Calculated total density of states for c-GST. The 
Fermi level is set to zero. 



FIG. 2: Calculated reflectivity of c-GST. Previous experimen¬ 
tal result from Ref.— is also shown for a better comparison. 



FIG. 3: Calculated optical absorption coefficient along 
with experimental data from previous measurements (see 

Refs.iSi^a) 
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FIG. 4: Calculated thermopowers of c-GST as a function of 
doping concentration both for p-type and n-type at various 
temperatures. Solid horizontal line at 200 t'V/K represents 
the limitation of Seebeck coefficients for good TE materials. 


the energy up to 3.5 eV which covers the energy range of 
DVD (1.91 eV) and blu-ray (3.06 eV)^ technology . This 
is consistent with the measured optical reflectivity by 
Garcia-Garcia et al.'s^ [Fig. [2], though they used thin 
Aims of c-GST in their experiments. Fig. [3] depicts the 
calculated optical absorption coefficients together with 
that from previous measurements. Our calculated ab- 
sorptioir coefficient compares well in magnitude with that 
from Lee et al^ and Pirovano et al^ for the higher en¬ 
ergy regime but disagrees at lower energy reflecting the 
band gap difference. One possible explanation is that 
the structure is more complex than the random struc¬ 
ture that we assumed, or that larger length scale effects 
such as local near ferroelectric polarizations play a role in 
the samples. Additionally, the dimension of the crystals 
may also play a role since thin films of GST were used 
in both of the experiments instead of the ideal bulk that 
we are considering. 


B. Thermoelectric properties 

In Fig. m we have plotted the calculated thermopower 
(S(T)) of c-GST at various temperatures ranging from 
100 K up to 500 K, for both electron-doped (n-type) 
and hole-doped (p-type) cases. It has been reported that 
it might be possible to dope GST n-type^i^, however, 
keeping in mind that most of the experiments deal with 
p-type GST, we focus on the thermoelectric properties 
of p-type c-GST. In fact, there is a rather good electron- 
hole symmetry in the S(T), so n-type if realized would be 
similar to p-type. We give the direction average of S(T), 
i.e. to eliminate artificial anisotropy due 

to the unit supercell. 

We find c-GST exhibits a Seebeck coefficients 200 - 
323 /iV/K when doped p-type for doping concentration 



FIG. 5: Temperature dependence of the thermopower of c- 
GST with two p-type doping concentrations and experimental 
results from previous work.— 

1.7x10^® - 6.1x10^® cm“® at room temperature. More 
importantly, at low temperature (T = 200 K), the See¬ 
beck coefficient exceeds 200 pV/K for a doping range 
of 1x10^® - 3.5x10^® cm“®. A strong bipolar effect was 
noted for T = 300 K and higher temperatures, however, it 
was not present at lower temperatures. This is connected 
with the band gap, which as noted may be too low in the 
present calculations. In order to assess the possibility of 
c-GST as a low temperature TE material, it is necessary 
to compare the calculated Seebeck coefficient with that 
in popular low-temperature TE materials. For example, 
Poudel and co-workers^ reported a bulk nano-composite 
Bia;Sb 2 -a:Te 3 with S(T) of 185 pV/K at 300 K. Chung 
et alA discovered CsB^Tee with a thermopower of 105 
pV/K at 300 K. Importantly, the high thermopowers are 
at high carrier concentrations where reasonable electric 
conductivity may be anticipated. 

We shall now focus on comparing our calculated re¬ 
sults with previous experiments on GST thin films. Kato 
and Tanaka— investigated thermoelectric properties of 
Ge 2 Sb 2 Te 5 thin films with sample dimensions ranging 
from 50 nm - 600 nm for a temperature range 25 °C to 
200 °C. For their c-GST samples, they found a pretty 
high values of S(T) from 280 pV/K to 380 pV/K. Lee et 
oL— also measured the temperature dependence of S(T) 
with 25 nm and 125 nm Ge 2 Sb 2 Te 5 thin films from 20 
°C to 300 °C and observed S(T) as high as 205 pV/K 
and 175 pV/K for 125 nm and 25 nm thick films, respec¬ 
tively. Both of the authors found a decreasing S(T) with 
increasing temperatures. In order to further compare our 
calculations with these measurements, we matched the 
calculated Seebeck coefficients at 300 K with that from 
Kato’s data for the 550 nm thick thin film and calculated 
the corresponding doping levels in c-GST. For these par¬ 
ticular doping levels, we then plotted the variation of 
S(T) with temperature together with the respective ex¬ 
perimental data, as shown in Fig. [5] As seen, a similar 
decreasing trend of S(T) with respect to temperatures 
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can be found in the calculation for both doping concen¬ 
trations. However, the agreement of the decreasing fea¬ 
ture is more obvious in the low doping level case which 
falls into the bipolar region and the deviation may stem 
from the underestimation of the band gap. This may also 
suggest the fact that the experimental samples are in the 
low carrier concentrations range where bipolar effect is 
strong and detrimental to the thermoelectric properties. 
Therefore a better thermoelectric performance might be 
feasible in higher doping concentration samples. 

IV. SUMMARY&CONCLUSIONS 

To summarize, we performed first-principles calcula¬ 
tions to investigate electronic, optical and thermoelectric 
properties of c-GST with fully random cation site oc¬ 
cupancy generated using the SQS scheme. Our results 
based on Boltzmann transport calculations predict sub¬ 
stantially high thermopower at high carrier concentration 
(~ 300 ^V/K at 7x10^® cm“®) at room temperature. 


We also find high thermopower at high carrier concen¬ 
trations going down to 200 K and below suggesting that 
with heavy hole doping c-GST may be an excellent low- 
temperature thermoelectric. These results suggest future 
experiments to better understand the low-temperature 
thermoelectric properties of c-GST. Specifically, it would 
be of interest to study the temperature range of 200 K 
- 300 K with high doping concentrations ranging from 
1x10^® cm“® up to 6x10^® cm“®. 
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